%  TEST CASE BY HAND

it      = 1; 
keff    = 1.0;% alt_keff;%1.039435608478710;%1.050277290361834;

%  data format
%    m1
%    g1 -->  D_1  SigR_1  nuSigF_1  xi_1  s_1->s_1 s_2->s_1 ...
%    g2 -->  D_2  SigR_2  nuSigF_2  xi_2  s_1->s_2 s_2->s_2 ...
%    ...
%    m2 ...

dat = [ 1.500   0.030   0.000  1.00   0.000  0.200  
        0.420   0.080   0.035  0.00   0.200  0.000
        2.000   0.040   0.000  0.00   0.000  0.000
        0.300   0.010   0.000  0.00   0.040  0.000 ];    
numg = 2;  
numm = 2;
xcm  = [   0   10  20 ];
xfm  = [     2    2   ];
ycm  = [   0   10  20 ];
yfm  = [     2   2    ];
mt   = [ 1 1
         1 1];       
src(:,:,1)  = [ 1 1
                1 1];
src(:,:,2)  = [ 0 0
                0 0];
BCL = 0;
BCR = 0;
BCB = 0;
BCT = 0;           

IBSL = [ ...
    0.0297
    0.0297
    0.0297 ]*0;IBSL = ones(xfm,1)*0;%ones(xfm,1);%CRL;
IBSR = [ ...
    0.0467
    0.0467
    0.0467 ]*1; IBSR = ones(xfm,1)*0;%0*CRR;%(end:-1:1);
IBSB = [ ...
    0.0297
    0.0408
    0.0467 ]*1; IBSB = ones(xfm,1);%0*CRB;%(end:-1:1);
IBST = [ ...
    0.0297
    0.0408
    0.0467]*1; IBST = ones(xfm,1)*0;%0*CRT;
